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Abstract Low-frequency solar and interplanetary radio bursts are generated at 
frequencies below the ionospheric plasma cutoff and must therefore be measured 
in space, with deployable antenna systems. The problem of measuring both the 
general direction and polarization of an electromagnetic source is commonly 
solved by iterative fitting methods such as linear regression that deal simultane- 
ously with both directional and polarization parameters. We have developed a 
scheme that separates the problem of deriving the source direction from that of 
determining the polarization, avoiding iteration in a multi-dimensional manifold. 
The crux of the method is to first determine the source direction independently 
of concerns as to its polarization. Once the source direction is known, its direct 
characterization in terms of Stokes vectors in a single iteration if desired, is 
relatively simple. This study applies the source-direction determination to radio 
signatures of flares received by STEREO. We studied two previously analyzed 
radio type III bursts and found that the results of the eigenvalue decomposition 
technique are consistent with those obtained previously bv lReiner et al.\ (Solar 
Phys. 259, 255. l2009|) . For the type III burst observed on 7 December 2007, the 
difference in travel times from the derived source location to STEREO A and B 
is the same as the difference in the onset times of the burst profiles measured 
by the two spacecraft. This is consistent with emission originating from a single, 
relatively compact source. For the second event of 29 January 2008, the relative 
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timing does not agree, suggesting emission from two sources separated by 0.1 
AU, or perhaps from an elongated region encompassing the apparent source 
locations. 

Keywords: Solar wind, Direction finding, Type III bursts, Solar radio 



1. Introduction 

Solar activity encompasses a wide range of phenomena that may vary in complex- 
ity. Solar flares and coronal mass ejections (CMEs) are among the best known 
and most spectacular examples of this activity. During flares and CMEs, energy 
stored in the coronal magnetic field is converted into kinetic energy and thermal 
energy which are observed as coronal mass movement, and used to heat the 
plasma adjacent to the reconnection region, and in the acceleration of particles 
that eventually will follow magnetic field lines in the interplanetary medium. 

Flares and CMEs emit radiation at a variety of frequencies, from X-rays at 
the coronal sources to metric and decamctric radio emission in the interplanetary 
medium. Each observed frequency range can then be associated with a specific 
mechanism of emission, which gives clues about the physical processes taking 
place in the solar atmosphere or interplanetary medium when the emission 
emanated. 

The radio emission observed at metric and decametric wavelengths consists 
of a long series of "radio bursts" that may last from a period of hours to days 



(Bas tian, Benz, and Gary, 1998 ). The physical processes at work, including the 
changing ambient conditions, are signified by temporal evolution in the frequency 
spectrum of the emission. Historically, solar radio bursts have been classified 
into types, based on the morphology of features that appear in plots of elec- 
tromagnetic power in a frequency-time plane. Type II radio bursts are generally 
observed during major flares and/or CMEs associated with filament eruptions, at 
meter wavelengths for periods from 5 to 30 min. These bursts show a distinctive 
frequency drift in time, from higher frequencies to lower frequencies. Type III 
bursts are the most common, believed to be the result of flare-energized electrons 
that produce "plasma emission" in the meter and decimeter range with lifetimes 
of approximately 10 s. In general, type III events appear in groups lasting several 
minutes, which indicates that each observed event is a collection of many small- 
scale type III bursts. In this paper we are interested in interplanetary radio 
emissions from type III bursts, particularly the direction of the radio source, but 
the results can also be applied to type II bursts. 

Radio direction finding (DF) is a powerful technique that allows us to de- 
termine the direction of a radio source, using observations made with simple, 
deployable electric antennas. There are several methods to do this, depending 
on how the radio emission is received and recorded. Some methods require distin- 
guishing signal modulation in receivers on spinning spacecraft ( |Fainberg, Evans, and Stone, 1972] ) .| 
Others are based on auto- and cross-correlation analyses of data from three-axis- 
stabilized spacecraft, such as on NASA's Solar Terrestrial Relations Observatory 
(STEREO) ( |Kaiser et al, 2008| [Bougeret et al, 2008D , Cassini, and future mis- 



sions like Solar Probe Plus and Solar Orbiter. Analysis techniques for these 
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data sets can be divided into two groups: iterative fitting ( |Ladreiter et al.~1 995 



Vogl et al., 2004| Reiner et al., 2009) and noniterative methods ( SantoHk, Parrot, and Lefeuvre, 2003[ [ 



Cecconi and Zarka, 2005 |Hess, 20T0| ). To fit observational data, the methods 
in the first group use theoretical models of the auto- and cross-correlations to 



find the direction and polarization of the radio emission (Ladreiter et al., f995 



Cecconi and Zarka, 2005 1 . The method used bv lReiner et aZ.I(|2009|) and Vogl et aH 



( 2004J) uses a least-squares minimization fitting of the data, while the method 
described in lLadreiter et aLl(|1995p iterates the solution obtained using a singular 
value decomposition (SVD) process. The ILadreiter et al.\ (|f 995[) method is a 
more sophisticated iterative fitting method that allows one to estimate the factor 
of error of each solution found by the iteration process. 

Noniterative methods encompass different techniques to solve the direction 
finding problem. For most polarizations, it is possible to find the solution by 
inverting the particular case of the general solution of the auto- and cross- 
correlations, as is shown in Cecconi and Zarka (2005). This method gives an effi- 
cient and fast solution to the problem; however, its downside is that it constrains 
the polarization of the radio signal, which is generally unknown. SantoHk, Parrot, and Lefe uvre^ 
(|2003p implemented a method that locates the direction of emission of electro- 
magnetic waves, by finding the singular values of the spectral matrix, a 3 x 6 
array composed of the values of the auto- and cross-correlations, noting that the 
values of the cross-correlations are generally complex. The relation between the 
singular values obtained using an SVD algorithm gives the direction of the radio 
source and certain properties of the polarization. 

All the foregoing examples address the problem of simultaneous determination 
of both the source direction and its polarization, the latter in terms of Stokes vec- 
tors. In this paper we develop a simplified method to derive the general direction 
of the radio source distribution independently from its polarization. This method 
is similar in form to the SVD method proposed by |SantoHk, Parrot, and Lefeuvre 



(20033, but it addresses only the problem of finding the location of the radio 
source, using an eigenvalue decomposition method. We apply the method to two 
type III bursts analyzed by IReiner et al. \ (|2009[) and compare our results with 
theirs. 



2. The Eigenvalue Decomposition Method 

We will suppose that the wave electric field E in the vicinity of the spacecraft 
can be expressed at a specific time by 

E r (t) = Rc{Ee Mt }, (f) 

with the harmonic wave frequency, w, assumed constant, and with E generally 
complex. We first consider the case in which E is constant, noting the following: 

1. E r (i) moves in a plane, V , generally tracing out an ellipse. 

2. If E r (t) is the result of a single point source, then the direction of the source 
from the spacecraft must be along the line, 1Z, perpendicular to V , if V is 
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uniquely determined. This could be either of both directions from the space- 
craft in which 1Z extends. If the motion of E r (£) is confined to some line, C, 
then there is no unique V, in which case the source can be anywhere in the 
plane perpendicular to C in which the spacecraft sits at any moment. 

To prove that E r (i) moves in some plane, simply observe from Equation (fTJ) 
that 

E r (f ) = Rc{E} cos u:t - Im{E} sin ut , (2) 

i.e., E r (£) is a linear combination of two vectors Rc{E} and Im{E}. If Rc{E} and 
Im{E} arc linearly independent, then they uniquely define a plane, V , and E r (i) 
must be contained in that plane. Otherwise, if Re{E} and Im{E} are parallel, 
E r (£) is confined to a linear motion, which is contained in an infinite number of 
planes. Either way, motion in at least one plane is secured, but if Rc{E} and 
Im{E} are linearly independent, then the plane, V , in which E r moves is unique, 
whereas if they are parallel, then V is indefinite, and the direction to the source 
accordingly ambiguous. 

We defer the proof that the motion of E 1 in planeV is elliptical, as a deter- 
mination of source direction does not depend on it. 

The second point above is a result of well-established classical radiation 
theory, described at length in most texts. It applies to any source, not just 
a monochromatic one such as Equation JTJ) expresses when E is constant. 

When Re{E} and Im{E} are linearly independent, the source direction in the 
case of an electric field, E, known to be from a single source can be expressed 
by the unit vector, R, in the direction of 

R = Rc{E} x Im{E}, (3) 

meaning that the source is either in the direction of R or — R. Otherwise (i.e., 
Rc{E} and Im{E} are parallel), R is null, in which case the source can be 
anywhere in the plane perpendicular to either. (Besides the above, there is only 
the extreme case of a null source, in which both Re{E} and Im{E} are null. In 
that case, of course, the source can be anywhere at all.) 

2.1. The Case of a Distributed Source 

In the case of a completely coherent distributed source, meaning that the resul- 
tant E from all of the source components is temporally invariant throughout the 
measurement, E r (t) would remain confined to move in a plane, as proved above. 
In practice, far-separated source components tend not to be fully coherent. Con- 
sequently, Equation ([T|) can still express the motion of E r (i), but the direction, 
amplitude, and phase of E are now all subject to some variation, such that E r (t) 
is no longer strictly confined to a plane. For this problem, it is useful to regard R 
as the direction through the center of the source distribution by some definition. 
If the bandwidth of the radiation received by the spacecraft is Aw centered on 
w, then the plane in which E r (t) moves will itself tend to nutate with respect to 
R on a time scale of 1/Aw. Hence, E r ■ R will no longer be entirely null. 
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If the source is relatively compact, then E r • R will always be relatively small. 
Wc therefore now define a representative "source direction," R, for a distributed 
source to be cither of the two opposite directions from the spacecraft for which 
((E r • R) 2 } is minimal (or infinitely many directions for the various cases in 
which a unique minimum is nonexistent). Here, the angular brackets signify an 
average over the time interval of the measurement. The problem, then, reduces 
to the determination of the line TZ on which R lies when ((E 1 • R) 2 } is uniquely 
minimal or the plane to which R is confined, if any, when ((E r • R) 2 ) is not 
uniquely minimal. 

The relative magnitude of the minimum of ((E r • R)"} with respect to the 
total power, ((E r ) 2 ), can be regarded as the square of a characteristic radius, 
r c , of the source distribution. This is explained in Appendix |A"1 

2.2. Observational Data 

In general, spacecraft have a reference frame, S, that can be expressed by a set 
of orthonormal vectors, 

S = {ei,e 2 ,e 3 }, (4) 

pointing in the x-, y-, and z-directions, respectively. Hence, E can be expressed 
in terms of complex components, E a , thus: 

E = Y,E a e a , (5) 

a 

with 

a = {1,2, 3}. (6) 

For STEREO, the products used to determine the position of the radio source 
are complex statistical correlations between electric-vector components for some 
frequency, u, comprising a 3 x 3 matrix, C, whose elements are 

C a f) = (E a Ep), (7) 

with a, f3 € {1,2,3}, in some bandwidth, Aco, about w, where E is defined by 
Equation ([T]) and the angular brackets indicate an average over a time interval, 
At. 

2.2.1. Interpretation of the Correlations 

It is fairly straightforward to determine ((E 1 • R) 2 ) in any direction from C in 
any reference frame. We first define the real correlation, C r , by 

Clp = (Elm^t)), (8) 

where the angular brackets indicate averaging over time, and show that this can 
be determined from C: 

Clp = (Re{E a e iut }Re{E^e iut }) 

= \({E a e 1 ^ + EZe-^XEpe*** + ^e"^*)). (9) 
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At this point, we will apply the distributive property to the product on the right 
side of Equation noting that this will have terms containing exp(mwt) with 
n G {0, ±2}. We now assume that the spectrum of E outside of a bandwidth 
Auj < u>/2 is either null or negligible. Thus, we can drop all terms for which 
n G {±2}, since the averages of these will be null (or similarly negligible). This 
leaves us with only the cross terms in which an occurrence of exp(zu;£) is offset 
by its inverse, exp(— iut), hence, 

C a/3 = j{ E ocEp +EpE* a ) 

= hi C {(E a E;)} 

= hic{c a p}. (io) 

Stated slightly more succinctly, 

C r = iRc{C}. (11) 

Thus, the correlations between the real electric components averaged over time 
are simply half of the real parts of the complex correlations. 

Note that the diagonal components of C r are the mean squares of the electric 
components, ((E 1 • e Q ) 2 ), in the frame of the spacecraft. Now, given the real 
matrix, C r , in any orthonormal coordinate system, S, it is straightforward to 
determine its analog, C r ', in any other orthonormal frame, 



S' = {&[, e' 2 , e> 3 }, (12) 

from its diagonal components ((E r -e^) 2 }, in S' . Hence, knowing C r in any single 
frame, S, allows us to determine C r/ in any other frame, 5', thus to take from 
its diagonal components ((E 1 • e' a ) 2 }, in S'. 

To see that this is true, let U represent the orthonormal transformation from 
S to S' . To be particular, we express U by the real, orthonormal matrix whose 
elements are U a p, with a, /3 G {1,2,3}, in that the components, E' a = E • e' a , 
of E in S' are 

ei' = y. u ^ e I- ( 13 ) 

Extending Equation (8) to the transformed reference frame, it now follows that 



^2 ^c^E U PvEl 

/Z UafiiE^E^Ujp 

fit' 

Y.U^Cl^ (14) 
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where Ujg = Up v defines the transpose of U, and because U is real and orthog- 
onal, U T is the inverse of U. Hence, 

C r ' = UC r U _1 . (15) 

Note that, because multiplication is commutative, 

(E r a Ej) = (EjEl), (16) 

and it therefore follows that 

C a /3 = Cp a , (17) 

i.e., C is symmetric. Thus, there exists an orthonormal rotation, U, that di- 
agonalizcs C 1 . This is simply the matrix whose columns are the eigenvectors of 
C r . 

Since C r is real and symmetric, its eigenvalues must be real. Indeed, since the 
diagonal components of E V J are correlations of components of E r with themselves 
for whatever U, the eigenvalues must be non-negative. These eigenvalues will 
generally have a greatest and a least value. Broadly speaking, a three-dimensional 
map of ((E r -f) 2 ) as a function of direction, f , traces a three-dimensional ellipsoid 
whose minor axis is the least eigenvalue of C r in length, in the direction of the 
corresponding eigenvector; whose major axis is the greatest eigenvalue of C r in 
length, in the direction of the eigenvector corresponding to that (perpendicular 
to the first); and to whose intermediate axis (perpendicular to the other two) 
is attached an eigenvalue somewhere in between the least and the greatest. The 
source direction, R, is taken to be along the minor axis of the foregoing ellipsoid. 

There are special cases leading to degeneracies in which any two or all three 
eigenvalues can be equal, or in which one or more can be null (all three for a 
null source). In those cases, the source direction becomes less definite in certain 
respects similar to those previously noted in the special case of a point source. 
In the case of two equal eigenvalues less than 1/3, we can only suppose that R 
lies in the plane spanned by the two corresponding eigenvectors. In the extreme 
case of a null source or, e.g., a completely isotropic distribution of randomly 
polarized noise, all three eigenvectors are equal, and the source direction is then 
completely indefinite. When there is a single least eigenvalue, the eigenvector 
that corresponds to it is the one to be identified with R, the direction toward 
(or away from) the center of the source distribution. 

Thus the problem of finding R reduces to the eigenvalue problem of diag- 
onalizing the real matrix Re{C}/2, selecting the least of its eigenvalues and 
identifying the eigenvector that corresponds to that eigenvalue. The crux of the 
eigenvalue problem is solving for the roots of a secular equation that is generally 
cubic. All three roots are ensured to be real and positive. The cubic equation 
can be solved in any number of ways. 

Once the source direction is determined from a minimal eigenvector, one can 
straightforwardly determine the polarization of the source in terms of Stokes 
vectors from the directions of the maximal and intermediate eigenvalues (or 
multiply maximal in the case of one or more equalities) and the eigenvectors 
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2007 Dec 07 Time 

Figure 1. STEREO A and B dynamic spectra of the 7 December 2007 type III burst from 
04:00 UT to 06:00 UT. The plotted frequencies range from 100 kHz to 16.025 MHz. The color 
bars represent the intensity of the radio emission measured in decibels. 

attached to them. This is the subject of a continuation of this study. The relative 
value of the minimal eigenvector similarly allows us to model the effective angle 
subtended by the source distribution. 

3. The STEREO Radio Observations 

For application of the formalism to STEREO/WAVES observations, we devel- 
oped a set of codes that allow us to determine and plot the positions of type II 
and type III radio sources. In this section we present and compare the results of 
the method applied to two well-studied type III bursts bv lRciner et aLI(|2009[) . 

3.1. 7 December 2007 Type III Event 

A solar flare classified as a GOES B1.4 was observed on 7 December 2007, hosted 
by the active region 10977 located close to the disk center (S05W06). The flare 
started at 04:36 UT, reaching maximum at 04:39 UT and ending at 05:01 UT. 
This flare was associated with a type III burst registered by space- and ground- 
based radio observatories, including both STEREO and Wind spacecraft, with 
a starting time of about 04:36 UT (Figure [1]). 

During the period of observation STEREO A was at about 20.8° of separation, 
in the Earth direction of rotation and at a distance from the Sun of 0.967 AU, 
while STEREO B had a separation of about 21.6° in the opposite direction and 
a distance of 1.027 AU from the Sun. 
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Figure 2. Results of the direction finding analysis at STEREO A and B at 325, 425, 625, 875 
and 1075 kHz for the 7 December 2007 solar type III burst. The top panels show the derived 
radio emission power, the middle panels show source azimuths with respect to the line of 
sight between the Sun and the respective spacecraft, and the bottom panels show the relative 
elevation to the ecliptic plane. The horizontal black dashed lines show the power threshold 
used to discriminate the measurements to be plotted. 
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As mentioned in the previous section, we characterize source directions with 
respect to either spacecraft in terms of an azimuth and an elevationQ In this 
context, the azimuth is dehncd as the angle between the STEREO-Sun line and 
the source direction. Similarly, the elevation is defined as the angle between the 
ecliptic plane and the source. Figure [2J shows the direction finding calculations 
for various frequencies for each of the STEREO spacecraft and also shows that 
the source azimuth decreases with frequency from the perspective of STEREO 
A. This parallax suggests sources distributed approximately along a line of sight 
from STEREO B with the higher frequency sources closer to the Sun so that 
they have a greater azimuth with respect to STEREO A. Also, a clear time delay 
is observed in the top panels in Figure [2J the implications of this delay will be 
discussed later. This figure shows the practicality of the eigenvalue decomposition 
method described in Section [2J 

The 425 kHz signature was chosen to compare our calculations to lReiner et al.\ 
( 200SJ) , who use an iterative method to determine the locations of the radio 
sources, based on theoretical relations of the antenna power. lCecconi et q/.l(|2008p [ 
describe the theoretically deduced auto- and cross-correlations, and the angles 
indicating the direction of arrival of the radio emission. 

Table[T]shows the results obtained by the two methods. The average values of 
azimuth calculated around the maximum signal for the 7 December 2007 event 
are about -1.7° for STEREO A and -0.2° for STEREO B; the elevation angles 
are —5.2° and —0.7°, respectively. These results are close to the values reported 
in lReiner etaJ\ ([20051) (-4.5° of azimuth in the reference frame of STEREO A 
and +3.0° in the frame of STEREO B). Although there is a small discrepancy 
between our azimuths and those of lReiner et aZ.I (|2009|) . the results appear to be 
consistent with a Wind /WAVES observation that shows a source azimuth from 
a near-Earth perspective of —0.5°, also reported bv lReiner et al.\ (|2009p . 

Table 1. Comparative table of the values obtained for 
the 7 December 2007 type III event, using a least-squares 
fitting (LSF) algorithm and the eigenvalue decomposition 
(EVD) algorithm. 





STEREO A 


STEREO B 


Wind 




LSF 


EVD 


LSF 


EVD 


LSF 


Azimuth 


-4.5° 


-1.7° 


+3.0° 


-0.2° 


-0.5° 


Elevation 


« 0°„ 


-5.2° 


« 0°, 


-0.7° 





* Reported as "close to zero." 



The positions of the 425 kHz radio sources relative to the Sun for each tempo- 
ral observation are plotted in Figure [3l The coordinates of the center of the solar 
disk are (0°,0°). The closed curves in Figure [3] represent circles on the celestial 
sphere centered upon respective source centroids and whose radii, r c , are the rms 



1 The elevation is denned as the angle measured between the ecliptic plane and the source 
direction. [Reiner et a/.M2U09l makes use of co-latitude instead of elevation, which has a simple 
relation that is defined as elevation = 90° — colatitude 
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Figure 3. Azimuth-elevation plot of the derived positions at 425 kHz for the 7 December 
2007 solar type III burst. The concentration of points close to the origin represents the gen- 
eral location of the type III radio source. Data points signifying directions far from the Sun 
characterize pre- or post-flare instrumental and background noise. Also plotted for each such 
direction is a circle on the celestial sphere centered upon it whose radius is the characteristic 
radius, r c , of the respective location on which it is centered (see Appendixf^. Note that circles 
on the celestial sphere are warped when projected into an azimuth-elevation format. 



radii of a distribution of randomly phased sources so centered and consistent with 
the radio signatures (see Appendix [3} . Note that circles on the celestial sphere 
are generally warped when projected onto a planar azimuth-elevation format. 
(Indeed, circles that separate the north and south poles project into curves that 
are not even closed in an azimuth-elevation mapping.) This figure shows the 
source circles and their centroids for an extended period around the type III 
burst (04:30 UT to 05:00 UT), which include those for the background signal, 
which are centered in the southern hemisphere, much weaker, appear to be highly 
diffuse, and not necessarily the same for STEREO A as for STEREO B. 

The location of the source in the interplanetary medium was estimated by 
parallax, i.e., "triangulation," based on lLiu et al.\ (|2010p using the mean values 
of linc-of-sight azimuths and elevations plotted in Figure [2) A brief summary of 
the trigonometry for this application is given in Appendix [Bl 

The results of the parallax determination for the 7 December 2007 event are 
shown in Figure 0J for the specific case of the 425 kHz frequency the sources 
are about 28° east and 3° south, at a distance of 0.04 AU (triangle), which 
is not consistent with the result found bv IReiner et al.\ (|2009[) . of 2° east and 
pa 0° of elevation, at a heliocentric distance of 0.2 AU (square). This can eas- 
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ily be explained by the fact that a small change in the angular configuration 
may lead to large discrepancies in the parallax determination. The parallax 
algorithm works by finding the location of the intersection of the two vec- 
tors defined by the average positions. The fact that these locations are dif- 
ferent, even by a small amount, means that the intersection point can be in 
a considerably different location, as demonstrated above. The algorithm was 
tested using the results from lReiner et a l. (2009), obtaining values in the general 
neighborhood of those reported. Using the electron density model described 
in |Leblanc, Dulk, and Bougeret| (|1998p with an electron density at 1 AU of 
7.2 cm -3 , we calculated the distance where the radio emission could be produced 
for a range of frequencies. We found that, for this model, a wave of frequency 
425 kHz is emitted at a distance of about 0.058 AU, which is within the limits of 
the values obtained by parallax, therefore showing the plausibility of our results, 
assuming that the radio emission was generated at the fundamental plasma 
frequency^] 

Figure [4] also shows the canonical Parker spiral, representing the likely prop- 
agation path of the emitting electron beam generated by the flare. However, 
remember that electron beams have a physical size and are likely to be in- 
homogeneous. This, along with changing conditions in the media into which 
they impinge, affects the location of the centroid of emission, which is what 
this direction finding technique and others want to determine. The frequency 
drift of the radio sources is seen in Figure 0] high-frequency sources are located 
closer to the Sun, while low-frequency sources are further in interplanetary space. 
However, bear in mind that computed positions at higher frequencies may be 
unreliable due to uncertainties in the antenna pattern (see the discussion at the 
end of this section). 

We now propose to examine the timing of the radio profiles at the two space- 
craft as a prospective control resource for the locations determined by parallax. 
We do this by first computing the distance to each spacecraft from the extrapo- 
lated locations. Then, times of flight for each spacecraft are computed assuming 
that the radio emission travels in a straight line from the source centroid to the 
spacecraft at a constant velocity, the speed of light. The difference between these 
two times is compared with the time shift between radio flux profiles at the two 
spacecraft (sec top panels of Figure [5]) to determine whether they are consistent 
with the source locations extrapolated by parallax. This "time-of-flight analysis" 
assumes that the onsets of the signals at the two spacecraft are the signature 
of radio emission simultaneously emitted from a single compact source. The 
limitations of this "time-of-flight analysis" are in the temporal resolution of 
the measurements and errors inherent in the assumptions of compactness and 
simultaneitjjfl 



2 It should be noted that it is still unknown whether type III radio emission is produce at the 
fundamental or harmonic plasma frequency, or both. 

3 The real radio sources appear to be spatially extended, and emission from spatially extended 
sources can be directional, emanating toward one spacecraft at one time and the other at 
another time. 
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Figure 4. Radio source locations are plotted in a Sun-centered reference frame with Earth 
at zero azimuth. In this format, source azimuths are the Earth-Sun-source angle projected 
onto the ecliptic plane. The horizontal axis in the top frame represents the ecliptic with the 
Sun at the origin. Bottom: Locations of the radio sources determined by parallax projected 
perpendicularly onto the ecliptic plane. The 425 kHz source location is indicated by a triangle 
data point. The location of the 425 kHz source reported by Reiner et al. ( 2009 ) is indicated by 
a square data point. The Parker spirals are plotted in gray and calculated using the formula 
(p = 0o — (©0 /V Bvi )r, where r, is the radial distance to the Sun, <j>0 is an arbitrary angle, 00 is 
the rotational velocity of the Sun (2 km s _1 ), and V sw is the solar wind velocity (400 km s~ 1 ). 
The cross symbols represent positions that are unreliable due to uncertainties in the antenna 
pattern at high frequencies (sec discussion at the end of Section 3.1) 
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We find that the time shift computed from the direction finding results ranges 
from «0.9 min at the lowest frequencies decreasing to «0.15 min at the highest, 
while the observed delay between the onsets of the emission received by STEREO 
A and B range from 1.3 to 0.65 min, respectively. For a comparison to other work, 
we consider the 425 kHz radio source. For this frequency, we find a time delay 
of 0.4 min for both our location and that of lReiner et al.\ (|2009p . although the 
former is significantly more distant than the latter. The time delay between the 
radio signatures is about 0.65 min. We found that the direction finding and 
timc-of-flight analysis results are consistent within the errors inherent to both 
techniques. 

For frequencies higher than «1000 kHz the radio signatures are unreliable, 
the parallax algorithm sometimes giving obviously unrealistic locations. The an- 
tenna pattern is more complicated at high frequencies, as the radio wavelengths 
approach the quarter- wave resonance of the antenna system ( |Bale et al, 20 08 ) , 
making the antenna calibration more difficult. Ongoing efforts to improve the 
calibration of the STEREO /WAVES instrument based on a better understanding 
of the effective antenna lengths and directionality will improve the results of 
the direction finding techniques in the future at high frequencies. The forego- 
ing uncertainties do not necessarily imply that the parallax determinations are 
generally wrong, just that they are unreliable. 

3.2. 29 January 2008 Type III Event 

The other type III event studied bv lReiner et aLI ( |2009[) occurred on 29 January 
2008. This specific type III was associated with a GOES B1.2 class flare, hosted 
by the active region 10982, which was located close to the east limb of the 
Sun. The flare began at 17:28 UT, reaching maximum at 17:34 UT and ending 
at 17:43 UT (Figure [5]) • The maximum of the radio emission was observed at 
17:45 UT ( |Reiner et al, 2009| ) by both STEREO spacecraft. On this date both 
STEREO spacecraft had a separation between them of 45.2°. STEREO A was 
located 21.6° west, at a radial distance of 0.967 AU, while STEREO B was at 
23.5° to the east of the solar disk at a distance of 1.0015 AU. 

Figure [5] shows the radio flux and source directions relative to both space- 
craft. The radio flux was an order of magnitude more intense at the location 
of STEREO B than at STEREO A. There is also a greater spread of source 
directions for STEREO A. For STEREO A, we find a mean azimuth around the 
maximum of power at 425 kHz of —11.4° and a mean elevation of about —6.3°. 
For STEREO B we find a mean azimuth of about —10.0° and a mean elevation 
of —0.6°. These results are in reasonable agreement with lReiner et aZTI (|2009p . A 
comparison is shown in Tabled 

Figure [7] plots the source locations, determined by parallax. We found that 
the 425 kHz source was located at about 73.0° east and at elevation of about 
15.0° south at a distance of approximately 0.21 AU (triangle). These results 
are almost in agreement with the results obtained by IReiner et al\ (|2009|) , of 
64.0° east and a heliocentric distance of 0.19 AU (square in Figure [7]), given 
that slight differences in the angles used in the parallax calculation can lead to 
substantial discrepancies. This gives a difference of about 9.0° in azimuth and 
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Figure 5. STEREO A and B dynamic spectra of the 29 January 2008 type III burst from 
04:00 UT to 06:00 UT. The plotted frequencies range from 100 kHz to 16.025 MHz. The color 
bars represent the intensity of the radio emission measured in decibels. 



Table 2. Comparative table of the values obtained for the 29 
January 2008 type III event, using a LSF algorithm and the 
EVD algorithm. 





STEREO A 


STEREO B 


Wind 




LSF 


EVD 


LSF 


EVD 


LSF 


Azimuth 


-11.0° 


-11.4° 


-8.5° 


-10.0° 


-11.5° 


Elevation 


«0°„ 


-6.3° 


0°, 


-0.6° 





* Reported as "close to zero." 



0.02 AU in radial distance between our results and the results of lReiner et all 
In their work, IReiner et al.\ (|2009l) do not give precise information about the 
elevation angle obtained by parallax, although they describe that the location 
of the source is close to the ecliptic plane. 

The time-of-flight analysis for this event using the direction-finding calculated 
data prescribes time delays from wl.7 min at the lowest frequencies decreasing to 
0.45 min at the highest. The radio profiles themselves show time delays from «1.4 
min for the lowest frequencies decreasing to wO.O at the highest (see Figure |BJ 
top panels). For a frequency of 425 kHz, we expect a time delay of 1.0 min for 
our location, compared with 1.1 min for that of lReiner et aLI(|2009p . The time 
delay in the radio flux is about 0.65 min, about the temporal sampling interval 
of the STEREO/WAVES instrument. The results of the time-of-flight analysis 
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Figure 6. Results of the direction finding analysis at STEREO A and B at 425, 625, 775, 
925 and 1075 kHz for the 29 January 2008 solar type III burst. The top panels show the 
derived radio emission power, the middle panels show the derived azimuthal angles of the 
radio source and the bottom panels show the relative elevation to the ecliptic plane. The red 
dashed lines show the enhanced STEREO A power for a better visualization of the signal time 
delay between the two spacecraft. The horizontal black dashed lines show the power threshold 
used to plot only the azimuth and elevation associated with the type III burst. 
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Figure 7. Location of radio sources in the interplanetary space associated with the type III 
radio burst of 29 January 2008, for all operational frequencies. The location we derive for the 
425 kHz source is marked by a triangle. The location of the same reported by IReiner et al.\ 
( 2009) is marked by a square. The Parker spiral is plotted in gray and calculated assuming 
a constant velocity of 400 km s _1 . Positions at high frequencies are expected to be incorrect 
(see Section I3.1II . However, they are not obviously implausible, as seen for the event of 2007 
December 07. 



for the 29 January 2008 type III radio burst confirm the results obtained by the 
direction-finding algorithm, showing that the locations obtained are plausible. 



4. Conclusions 

The radio direction finding technique is a powerful tool in the study of radio 
emissions in the interplanetary medium. Several different techniques have been 
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developed during the years to study type II and type III solar radio bursts 
(IFainberg, Evans, and Stonel IReiner et all 115551 IMP ILadreiter etaB 

i-3===!T5^-^— ^ s-= = 1 1 i : :— = : 



"1995l|Vogl et al.[ 120041 |Santohk, Parrot, and Lefeuvre[ [20031 ICecconi and Zarkal 



2005; Hess. 2010). For three-axis-stabilized spacecraft, such as STEREO, some 
of these techniques have been applied with some success, leading to new insight 
into hcliosphcric radio emission during flares and CMEs. 

We have developed a simple eigenvalue decomposition method to find the 
locations of radio sources based on observations from STEREO/WAVES. This 
method can be easily adapted to observations from other missions, such as Solar 
Probe Plus and Solar Orbiter. We have applied the method to two observational 
data sets already analyzed by IReiner et al.\ (|2009p with similar results. For the 
7 December 2007 type III event, the results of the algorithm are closer to the 
results obtained by the Wind spacecraft, with a deviation of w 3° from the 
results reported by IReiner et al. \ (|2009l) . A slightly smaller deviation was found 
for the 29 January 2008 type III event. The statistical uncertainty in the source 
directions determined by both EVD and LSF techniques appears to be about 
1.0° dCecconi et al, 2008[[Hess, 2010| ). The difference between the EVD and LSF 
source locations tends to be in this neighborhood, with significant exceptions. 

Based on ILiu et ai\ (|2010p , we used parallax to extrapolate the locations of 
the radio sources for both type III events studied by IReiner et al.\ (|2009|) , for all 
operative frequencies of the HFR1 receiver of the WAVES instrument. When the 
parallax is small, i.e., the heliocentric angle between the spacecraft is narrow, 
small errors of the order of a few degrees in the angular position of the source 
introduce large uncertainties in the distance of the source from either spacecraft, 
and hence in its location with respect to the Sun. Significant angular uncertain- 
ties are inherent in spatially extended source distribution, since these span wide 
angles that depend upon perspective. Modeling of radio signatures resulting from 
realistic source distributions conditions could give us insight into discrepancies 
between the simple models based upon point sources and the reality. 

A time-of-flight analysis was performed to check the direction finding solu- 
tions. In both cases, the trend in the time of flight is according to expectation for 
a type III burst emanating from a source region in the interplanetary medium 
above the east solar limb as viewed from Earth or either spacecraft. The burst 
emission arrives at the spacecraft at different times within ±2 min. We attribute 
the onset-time difference to the two spacecraft being at different distances from 
the source, hence resulting in different travel times. For both of the events we 
examined, the travel-time difference decreases with increasing frequency. For the 
geometry implied by parallax, this is consistent with high-frequency sources been 
closer to the Sun than low-frequency sources. For the 29 January 2008 event, 
the onset-time differences are somewhat less than those expected at the lowest 
and highest frequencies, assuming that the radio signals emanate simultaneously 
from the single point extrapolated by parallax and travel at the speed of light to 
the respective spacecraft. The difference between the extrapolated and measured 
time delays appears to be related to the extended spatial distribution of the 
sources, wherein the radiation received by one spacecraft can have emanated 
from a significantly different part of the source region than the other. Moreover, 
the lines of sight from STEREO A to the high-frequency sources pass closer 



SOLA: SoPH_DF_final.tex; 3 March 2013; 20:38; p. 18 



Determination of Electromagnetic Source Direction as an Eigenvalue Problem 



to the Sun than the radial distances extrapolated for the source centroids. This 
radiation is therefore subject to significant refraction, complicating the geometry 
of what part of the source distribution STEREO A sees and when it sees it. 



Appendix 

A. Modeling Radio Source Distributions 

In general, it is widely accepted that type III radio bursts are only weakly 
polarized, with some special cases presenting a degree of circular polarization 
up to 25% (ie.. lReiner et aOl2007p . Ideally, then, the real correlation ellipsoid 
of the type III bursts should be an oblate spheroid, with the two major axes 
equal. Whether or not this is exactly the case, we propose to characterize type 
III bursts with a size parameter that supposes that this is the case. For that 
purpose, it is useful to consider a model in which the radiation received by the 
receiver is an incoherent superposition of sources uniformly distributed over a 
cone of radius 0q, and to determine the shape of the spheroid for such a model. We 
begin by considering a reference frame with the receiver at the origin and where 
the z-axis is the axis of the cone. We first consider a source in the x-z-jAane and 
angle from the z-axis emitting radiation polarized in the x-z-planc of unit 
mean power at the receiver. In that case, 



{E 2 ) = cos 2 6>, 

{E 2 y ) = 0, and (18) 

{El) = sm 2 0. 

In the alternative case of an emitter in the same location but polarized in the 
y-z-plane, 

{El) = 0, 

{E 2 y ) = 1, and (19) 

{E 2 Z ) = 0. 

An unpolarized source at this location can be represented as an incoherent 
superposition of the above: 

{El) = cos 2 0, 

{E 2 ) = 1, and (20) 

{El) = sin 2 0. 

For this type of source, we note that 

{El)= S m 2 6 = r 2 , (21) 
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where 



(22) 



is the radius of the circle on the unit sphere that makes an angle 9 with the axis 
of the cone. Moreover, 



(E 2 ) ee {E 2 x + E 2 y +El) = 2. 



(23) 



For a distribution, D(£l) , of incoherent sources over the interior of a conical 
domain of solid angle, Q, whose half angle is 9q, 



{El) = 



r 2 D(fl)dQ, 



(24) 



<0 



while 



In other words, 



(E 2 ) = 2 / D(Sl)dSl. 

9<6o 



(E 2 ) 



(25) 



(26) 



where (r 2 ) is the mean square radius of the distribution from the conical axis of 
points on the unit sphere for angles less than 8q. Accordingly, we characterize 
the source distribution by the following characteristic radius from the centroid 
of the source distribution: 



:(E 2 )\ 1/2 

'(E 2 ) 



(27) 



B. Ecliptic-Plane Perspective of Source Locations 

The method we use to represent source locations in an ccliptic-plane-based ref- 
erence frame is taken from ILiu et ai\ (|2010[) who developed the trigonometry 
we use and explained its limitations. The method defines the source position in 
space in terms of the lines proceeding from the individual STEREO spacecraft 
in their respective directions to the supposedly single source. The intersection 
of these lines projected onto the ecliptic plane defines the ecliptic azimuth, a, of 
the source. If a a is the azimuthal angle between the Sun center and STEREO A 
or B {a G {A, B}), / representing the ratio between the heliographic distance to 
STEREO A (d A ) and STEREO B (rf B ), and 7 the known longitudinal separation 
between the two spacecraft, then the angle, /3a, between the source azimuth and 
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that of the spacecraft, for example, can be expressed by 

sin qa sinfas+T)— f sin a a sin qr t \ n j ^ n 

if cka > and «b < 



tan/3 



'A 



sin qa cos(qb+7)+/ cos qa sin as 

sin qa sin(QB —l) — f sin qa sin qb 

— sin qa cos(qb — 7)+/ cos qa sin qb 

sin qa sin (as +7) — / sin a a sin qb 

— sin qa cos(qb+7)+/ cos qa sin qb 



if a A > and a B > (28) 



if aA < and qb < 



Knowing the elevations, e CT , of the source with respect to the ecliptic, it is then 
straightforward to determine the distance of the source above the ecliptic plane. 
Errors in the determination of the azimuths and elevation can generally lead to 
different distances above the ecliptic plane insofar as the lines from the STEREO 
spacecraft toward the source fail to intersect. In this study, we take the elevation 
of the source to be the mean value of the observed elevation angles. 
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